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1. Introduction 

The Toda lattice has some very important applications in the theory of physics of 
nonlinear processes (see [1]). It is known the inverse scattering method allows one to 
investigate in detail the Cauchy problem for the Toda lattice in the different classes of 
initial data (see [1]-[15] and references therein). The last problem for the doubly-infinite 
Toda lattice 

a n — (b n+ i — b n ) , • = J^, a n = a n {t) > 0, , . 

b n = a 2 n - a 2 n _ ± , b n = b n (t), n = 0,±1,±2,.... 

with fast stabilized or steplike fast stabilized initial data is investigated in [l]-[9] (see 
also references therein) by the method of inverse scattering transform. However, this 
problem is not studied in the case of steplike initial data, where a n tend to zero as 
n — > +oo (or n — > — oo). 

In this paper we study the Cauchy problem for the system (1.1) with initial data 

a„(0) ->■ 0, b n (0) ->■ as n ->■ +oo, 

£|n|{M0)-l| + |& n (0)|}<oo. (1.2) 

n<0 

The solution is considered in the class 

IK(t)||c[ ,T] ->■ 0, ||6 n (*)||c[o,ri -> 0, as n ->■ +oo, 

ll<9(*)llc[o,n <0 ° > ^ L3 ) 

for arbitrary T > 0, where 

Q(*) = EW(M*)-i| + IM*)l) • 

ra<0 

Note, we cannot apply directly method given in [l]-[9] for the case inf a n > 0, 
because the Jost solution with the asymptotic behaviour on an +oo does not exist 
in our case. On the other hand, method of inverse problem is used (see [10]) in the 
case when Jacobi operator associated with (1.1) has the continuous spectrum [a, b] of 
multiplicity two. But this method cannot be used when the spectrum of the Jacobi 
operator has a continuous spectrum of multiplicity one and a discrete spectrum. 

The paper is organized as follows. In section 2 we formulate some auxiliary facts 
to the inverse scattering problem for the Jacobi operator associated with (1.1)-(1.2). In 
section 3 we describe the evolution of the scattering data of problem (1.1)-(1.2) . 

In the last section we prove existence of the solution of the problem (1.1)-(1.2) in 
class (1.3). 
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2. The scattering problem 

Consider Jacobi operator L generated in £ 2 (—00, 00) by the finite-difference operations 

( L y) n = a n-iy n ~i + b n y„ + 
in which the real coefficients a n > 0, b n satisfy the conditions 
a n — > 0, b n — > as n — > +00, 
\n\ {\a n - 1| + |6 n |} < 00 . 

n<0 

The interval [—2, 2] is the continuous spectrum of multiplicity one of operator L (see 
[16], [17]). Beyond the continuous spectrum, L can have a finite number of simple 
eigenvalues Hk(t), k = 1, ...,p. 

Let us formulate some auxiliary facts related to the inverse scattering problem for 
the equation 

(Ly) n = Xy n , n = 0, ±1,..., AgC (2.1) 

Many of these facts can be found in [16], [17]. 

Let P n (X) and Q n (X) be solutions of Eq. (2.1) with initial conditions 

P_x(A) = 0, P (A) = 1, 
Qo(A) = 0, Qi(A) = -. 

We denote by L semi-infinite Jacobi operator generated £ 2 [0, 00) by Eq. (2.1) as 
n > and the boundary condition y_ x = 0. This operator is completely continuous. 
Moreover, the spectral function p(A) of L represented [18] in the form 

P(A) = E 

A„<A 

where A n is the eigenvalue of L and (3 n is the norm of the eigenfunction corresponding 
to the A n . 

As is known from [18]- [19], the right Weyl function of the problem (2.1) has the 

form 

j-00 r — A 



or 

?-2 



n=l A « _ A 

where A„ ->■ as n ->■ 00. It follows from [12]-[13] that for A 7^ X k , k = 1, 2, Eq.(2.1) 
has Weyl solution 

^n(A)= Q n (A) + m(A)P n (A), (2.3) 
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"on the right semiaxis" (such that l^n (A)| < oo. ). 

n=0 

Suppose that T is the complex A-plane with cut along the interval [—2, 2]. In the 
plane T, consider the function 



/ \ A A 2 



choosing the regular branch of the radical so that y — 1 < for A > 2. We often 
omit the dependence of z(X) on A in what follows. Thus, in the formulas involving z 
and A, we always assume that z is as in the above equation. 

It is well known (see, for example, [20]) that Eq.(2.1) has a Jost solution represented 
in the form 

/ n (A) = a n z~ n (l + ]T A nm z'A . (2.4) 

V m<0 / 

The coefficients are given by 

i n = , b n = A n _i — A n+ i_i. (2.5) 

Oin+l 

Without restriction of generality we can suppose that X m G (—2, 2) for any 
m = 1, 2, As known [16], [17], for XedT, X 2 ^ 4, A ^ X m identity 



^(A) = a(A)/ n (A) + a(A)/ n (A) (2.6) 

holds, where the function a(A) can be regularly continued to T. Note also, a(A) can 
have a finite number of coinciding simple zeros outside the interval [—2, 2], because, 
these zeros constitute the discrete spectrum /i k , k — 1, of the operator L. 
Introduce reflection R(X) coefficient by the formula 



^ (A) " a(A) ' 

The function R(X is continuous for A G dT. Setting n = — 1 and n = in the identity 
(2.6) yields the expression 



1 f (X) + R(X)f (X) 

m(A) " -^i^xjTmr^x) (2 - 7) 

The norming constants M k (t) corresponding to the fJ,k(t) are given as 

oo 

M k ~ 2 = E fnM , k = i,..., P . 

n=—oo 

The set of quantities {R(X); fj, k ] M k , k — 1, is called the scattering data for 

the Jacobi operator L. The inverse scattering problem for L is to recover the coefficients 
a n , b n from the scattering data. 
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In solving the inverse problem, an important role is played by the Marchenko-type 
basic equation. Define 

^|>*V- + ^/ 8r ^-"A (2.8) 

where z k = z((i k ), k = 1, ...,p. 

Then A nm and a n involved in (2.4) satisfy the relations 

F 2n+m + A nm + A nk F 2n +rn+k = 0, m < n < , (2.9) 
a~ 2 = l + F 2n + £ A nk F 2n+k , n < . (2.10) 

k<0 

To reconstruct the operator L, we consider Eq.(2.8) which is constructed by the 
scattering data. We find A nm and a n from Eqs.(2.9) and (2.10), respectively, the firs one 
having a unique solution with respect to A nm .The coefficients a n and b n are definides 
for n < by (2.5). / n (A) for n < are defined by (2.4). From the formula (2.7) we 
obtain Weyl function m(A). The spectral measure dp(X) can be found by the formula 

dp(\ n ) = lim (A n - A)m(A), n = l,2, .... 

A — y An 

Using the approach in [12], [13], [19], we can reconstruct semi-infinite Jacobi operator L 
by its spectral measure dp(X). Therefore, we find a n , b n for n > 0. 



3. Evolution of the scattering data 

In this section we use the inverse scattering transform method to solve the problem 
(1.1)-(1.2). Let a n (t), b n (t) be a solution of the problem (1.1)-(1.2) satisfying (1.3). 
Consider the Jacobi operator L = L(t) associated with a n = a n (t), b n = b n (t). Jost and 
Weyl solutions, reflection coefficient, spectral measure now depend on the additional 
parameter t G [0, oo). 

Theorem 1. If the coefficients a n = a n (t), b n = b n (t) of Eq.(2.1) are solutions 
to problem (1.1)- (1.2) in the class (1.3), then the evolution of the scattering data is 
described by the formulas 

R(X,t) = J R(A,0)e( 2 " 1 - 2 )* (3.1) 

fi k (t) = n k (0), k = l,...,p (3.2) 

M k ~ 2 (t) = M fc - 2 (0)e( 2 * _1 -*0*, Zk = z ^ k ) , k = 1, ...,p. (3.3) 

Proof. System (1.1) is represented (see, for example [8], [13]) in the Lax form 

L = [L, A] = AL - LA, (3.4) 

where A — A(t) are Jacobi operator in £ 2 (— oo, oo): 

1 1 
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Since (3.4) implies that the family of operators L = L(t) are unitarily equivalent (see 
[5], [8]), the spectrum of L = L(t) does not depend on n and (3.2) is valid. 

Let / n (A, t)and ip n (\, £)respectively be the Jost and Weyl solutions of the Eq.(2.1) 
with the parameter t. Consider the indetity (2.6) with the parameter t. As follows from 
[8], [12] the function 4-0 n — {Aip) n is also a solution of the Eq.(2.1) with the parameter t. 
Appling the operator — A to (2.6), taking into account that the Jost solution f n (X, t) 
does not depend (see [8], on t asymptotically, we obtain 



- W)n = 0(A, t) + \ (Z- 1 - z) a(X, t)) f n (X, t) + 

(3-5) 

+ (d(A, t)-\{z- 1 -z)a{\ t))U\ t). 

On the other hand, we find 

± tPo _ (AP)a = 5^, | P ., _ (APU = _ _, . 

Since P n (X, t) and Q n (X, t) are linearly independent, the function j- t P n — (AP) n can 
be represented as 

jP n - (AP) n = A(X, t)P n + D(X, t)Q n . 



Setting n = — 1 and n = in the last relation, we find that 
A(X, t) = D(X, t) = a\. 

Therefore, 



f t P n - (AP) n = b -^P n + a\Q n . 
The same arguments are valid for solution Q n (X, t). Thus, we have the formula 

Now by the formula (2.3) with the parameter t we find that 
&n ~ (A^) n = (a?_M\ t) + ^) Q n + 

(3.6) 

+ (m(A, t) + ^m(A, t) - l) P n . 
Since L = L(t) is selfadjoint and bounded, j- t ip n — (Aip) n must satisfy the relation 

- {A^) n = 9(X, t)^ n . (3.7) 

Hence, we can represent the function ^ip n — (Aip) n as 

j t i> n - {A^) n = 9(X, t)Q n + 9(X, t)m(X, t)P n . (3.8) 
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Comparing this identity with (3.6), we have 

9(X, t) = a 2 _M\, + (3.9) 

Further, according to (2.6), (3.5), (3.7), 

0(A, t)a(X, t)Tn + 0{\ t)~aTX~J)f n = (d(\, t) + \ (z' 1 - z) a(X, t)) 

+ (d(A, t) - \ {z~ l - z) a{\, t)) f n . 

Since f n and /„ are linearly independent, so substituting (3.9) into the last indetity, we 
obtain 

d(A, t) + \ (z- 1 - z) a(X, t) = (a 2 _ 1 m(A, t) + a(A, t), 

a(A, t) - \ (z- 1 - z) a(A, t) = (a\m{\, t) + ^) a(A, t). 
From this relations, we get 

R(X, t) = (z- 1 -z)R(X, t), 
which imply (3.1). 

Now, let g n (^k, t) be a normalized eigenfunction of L. Since the eigenvalues 
fik, k — 1, p, of this operator are simple, we have 

j t 9n - (Ag) n = cg n . 

Taking the scalar products of g n with both sides of this equality in I 2 (— oo, <x>) and 
using HV'nll^^oo oo) = 1 an d A* = —A, we obtain c = 0. Therefore, 

j t 9n ~ (Ag) n = (3.10) 

On the other hand, if a normalized eigenfunction g n (^k, t) corresponds to the eigenvalue 
fik, then 

9n(lJ>k, t) = C k (t)f n (fJ, k , t). 

This implies that M%(t) = c 2 k (t). By virtue of (2.4), we find that 

j t 9n ~ (Ag) n ~ (c k (t) + Z -^^c k (t^j z~ k n 
as n — >■ — oo. Taking into account (3.10), we have 

c k {t) + Zj ^p^c k {t) = 

This equation implies the relation (3.3). 
The theorem is proved. 

Using Theorem 1, we obtain the following procedure for solving problem (1.1), (1.2) 
based on the inverse scattering transform method: Initial data (1.2) is given. Construct 
R(X, 0),//fc(0), Mfc(0), k — 1, ....,p . Calcullate R(X,t) ,fj, k (t) ,M k (t) using formulas (3.1)- 
(3.3). Construct a solution by solving the inverse problem by applying approach of the 
section 2 with R(X,0),fi k (0), M fc (0), k = 1, ....,p replaced by (3.1)-(3.3). 
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4. Solvability of the Cauchy problem for the Toda lattice 



In section 3, while constructing a solution to problem (1.1)-(1.2), we assumed that this 
solution exists in the class (1.3). Let us now investigate its existence. 

Theorem 2. The problem (1.1)- (1.2) has a unique solution in the class (1.3). 

Proof. Denote by B the Banach space of pairs of sequences y = (yi, n , V2,n)™ = _ 00 
for which the norm \\y\\ B = sup(|yi >n | + \y 2 , n \) + E \n\ (\yi, n \ + \V2, n \) is finite. Then 

n>0 n<0 

(see [21]) the set C ([0, T] ; B) of the continuous on an interval [0, T] with respect to 
the norm ||-|| B functions is the Banach space. 
Let as assume that 



x l,n = < 



a n (t) for n>0, 
a n (t) - 1 for n < 0, 



(4.1) 



x 2,n = Kif). 

Then system (1.1) is equivalent to the system 

X 1,T1 



% x \,n 1^2,71+1 ~ %2,n) + \ [ 1 — $n,\n\) ( x 2,n+l ~ x 2,n) i 



x 2 . 



X 



1,71 



X l,n-1 + 2 (l - <J„,| n |) (X hn - X 2 , n -l) , 



(4.2) 



where b n ^ m is the Kronecker symbol. 

Denote by F the operator generated the right-hand sides of system (4.2). Note, 
operator F is strongly continuously differentable in the space C ([0, T] ; B). 

Now passing to the integral equation in the standard manner, we find problem (4.2) 
with initial conditions 



Xl, n (0) 



a„(0) for n>0, 
a n (0) - 1 for n < 0, 



(4.3) 



ar 2) „(0) = 6 n (0). 
is equivalent to the equation 

x(t) = x(0) + f F(x(r))dT 
Jo 



(4.4) 



Applying the principe of compressed maps, we find that problem (4.4) on some interval 
[0, 6} has a unique solution x(t)with finite norm || cr(^) || o([o 6]- b) < 00 • Let us show that 
this solution can be extended to the entire positive semiaxsis. Assume the opposite. 
Then there exists a point t* G (0, oo)such that problem (4.2)-(4.3) has a solution 
x(t) = (xi tTl (t), X2, n (t))on the interval [0, t*) but lim = oo. It follows 

from [8], [13] problem (1.1)-(1.2) has a unique solution(a„(t), b n {t)) in C°° ([0, oo) ; M), 



Inverse scattering transform for the Toda lattice with steplike initial data 



9 



where M = £°° (—00, 00) © £°° (—00, 00). Hence, according to the (4.1) problem (4.2)- 
(4.3) has a unique solution x(t) = (xi >n (t), x 2 ,„(t)) satisfying 

\xi, n (t)\ + \x 2>n (t)\ < C 

for any t G [0, 00), where C does not depend on t . We integrate the system (4.2) over 
a interval [0, t}. Then, using the last inequality, after some simple transformations , we 
get 

||x(*)|| B <2||x(0)|| B + (4C + 4) f \\x(r)\\ B dr, < t < t* , 

Jo 

which, according to the Gronwall's inequality implies 

Mt)\\ B <2Mm B ^ c+4)t - 



Therefore, our assumption that ^ li m Q \\ x (t) \\b = 00 is n °t correct and problem (4.2)- 
(4.3) has a unique solution x(t) = (xi j7l (t), X2, n (t)) G C([0, T] ; 5) for any T > 0. 
Integrating the system (1.1) over a interval [0, t] and using (4.1), we obtain that problem 
(1.1)-(1.2) be uniquely solvable in the class (1.3). 
Thus, the theorem is proved. 
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